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Slow (logarithmic) relaxation from a highly excited state is 
studied in a Hamiltonian system with many degrees of free- 
dom. The relaxation time is shown to increase as the exponen- 
tial of the square root of the energy of excitation, in agreement 
with the Boltzmann-Jeans conjecture, while it is found to be 
inversely proportional to residual Kolmogorov-Sinai entropy, 
introduced in this Letter. The increase of the thermodynamic 
entropy through this relaxation process is found to be propor- 
tional to this quantity. 
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Study of the relaxation from highly excited states is im- 
portant not only in physico-chemical systems but also in 
biological systems. It has been reported that excitations 
of some protein molecules are maintained over anoma- 
lously long time spans Such behavior is relevant to 
enzymatic reactions, and in particular to their biological 
functions, because these molecules must maintain their 
excited state (brought about, say, by ATP) to be able to 
proceed from one process to the next. 

Typically, a Hamiltonian system with a sufhciently 
large number of degrees of freedom relaxes to equilib- 
rium. This relaxation process, however, is not neces- 
sarily fast. For instance, some Hamiltonian systems 
can exhibit logarithmically slow relaxation from excited 
states, called Boltzmann-Jeans conjecture (BJC), which 
was first noted by Boltzmann explored by Jeans [||, 
Landau and Teller 0| , and then has been studied in terms 
of nonlinear dynamics In the BJ conjecture, the re- 
laxation to equilibrium is required, but existence of chaos 
is not explicitly assumed. On the other hand, irreversible 
relaxation is often studied in relationship with chaos. In 
the present paper we study slow relaxation of a type in 
agreement with the BJC, by using a Hamiltonian dynam- 
ical system with a large number of degrees of freedom, 
and explore its possible relationship with chaos. 

We consider dynamics with the Hamiltonian 
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describing a set of pendula, where pi is the momentum 
of the i-th pendulum. The potential V is chosen so that 
each pair of pendula interacts only through their phase 
difference 
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Hence, the evolution equations for the momentum pi and 
the phase 9i are given by 
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This form of the interaction is chosen so that there is 
an attractive force tending to align the phases of all the 
pendula. The thermodynamic properties of this model in 
the iV — > oo, as well as the finite-size effect, are investi- 
gated in Ref. . Here we mainly discuss the case with 

= 10, which is sufficiently large to exhibit the thermo- 
dynamic properties, and the results to be discussed are 
valid for a larger system. The temperature of the system 
can be defined as T = {2K) /N, which is a monotonically 
increasing function of the total energy E. 

As studied in Refs. ||^-|^, each pendulum in this sys- 
tem exhibits small-amplitude vibration when the total 
energy is small, while, as the energy is increased, some 
pendula begin to display rotational behavior over many 
cycles. The relaxation of the rapidly rotating pendula 
into a vibrating assembly is rather slow, since their av- 
erage interaction with the assembly almost cancels out 
over the slow time scale of the assembly, due to the rapid 
rotation. Here we concentrate on such slow relaxation of 
the rotational mode to the vibrational modes. 

The slow relaxation from a highly excited state is in- 
vestigated systematically in the situation that the excited 
state is prepared by applying an instantaneous kick to a 
certain pendulum in the system. An example of the re- 
laxation process following this kick is depicted in Figj^, 
where it is seen that the kicked pendulum continues to 
rotate in isolation, maintaining a large energy. The re- 
laxation time tfi is defined as the interval of time from 
the kick to the point of which kinetic energy of the kicked 
pendulum first becomes smaller than K/N . 

The relaxation process depends on the energy i?o im- 
ported by the kick and also on the state of the other N—1 
pendula. Hereafter we call this assembly the 'bulk'. The 
bulk state is parameterized by the total energy E (or the 
temperature T) of the system before the kick | |ic| ]. 

The mean relaxation time (ii?) over an ensemble of 
initial conditions was computed by fixing E and varying 
Eq. As shown in Fig.0, the relaxation time increases 



exponentially with the kicked energy Eq | 1 1 | as 
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{tji) ocexp{a^/Eo), (4) 

for sufficiently large values of Eq. Here (•) represents an 
ensemble average and a is a constant. When the value of 
Eq is fixed, (tfl) decreases as a function of E, as is seen 
in the figure. 

The above exponential form of the relaxation agrees 
with that given by the BJC, which describes the relax- 
ation time for a system consisting of two parts with very 
different frequencies as an exponential of the ratio of the 
two frequencies. Assuming the separation of time scales 
of the kicked pendulum and the bulk pendula, this ex- 
ponential growth can be obtained through the Fourier 
analysis of the interaction term, as shown by Landau 
and Teller g] (see also g). 

In this argument, chaotic dynamics are not explicitly 
included, while chaos is often relevant to the irreversible 
relaxation to equilibrium. Actually, we have observed 
the exponential law (Q) in numerical studies clearly, as 
long as the Haniiltonian dynamics for the given energy E 
uniformly exhibits chaotic behavior without remnants of 
KAM tori. In this regime, (tn) decreases with the energy 
E, as shown in the inset of Fig.|[ On the other hand, 
in the regime of the small values of E, the relaxation 
process highly depends on the initial condition and it is 
difficult to define average relaxation time. The hindrance 
of relaxation is associated with the weakness of chaos, 
since the relaxation of the rotating pendulum stops when 
the dynamics of the bulk is trapped to the vicinity of 
remnants of tori. Now, we study a possible relationship 
between the relaxation time and the chaotic dynamics of 
the bulk. 

First let us consider the nature of the dynamics in the 
limit Eq — > oo. In this case, the kicked pendulum is com- 
pletely free from the bulk of vibrating pendula, whose 
collective motion is simply that determined by the Hamil- 
tonian of its {N — 1) pendula. At a finite value of Eq, the 
dynamics of the bulk deviate from those in the Eq —^ oo 
limit, due to the interaction with the kicked pendulum. 
In order to quantify such deviation, we choose bulk states 
with energies distributed around a given E, and compute 
the Lyapunov spectrum Ai(i?o) for various values of Eq. 
Note that, in the limit Eq — > oo, Xn-2 = correspond- 
ing to the free rotation of the kicked pendulum, while 
Aat-i = Xj\j ~ always holds due to the conservation of 
the total momentum and energy. 

Since the kinetic energy Kr of the kicked pendulum 
fluctuates and relaxes, Ai(i?o) is computed in the case 
that Kr satisfies \Eq — Kr\ < 6, with a small constant 
5. Since the relaxation of Kr is very slow, this condi- 
tion is satisfied over a time sufficiently long to compute 
the local exponents. In order to insure the convergence 
of the Lyapunov exponents, we take an ensemble aver- 
age of such local exponents (The existence of such 
time scale for the convergence is guaranteed by the slow 
relaxation, and indeed is confirmed numerically.) 



The values of the computed Lyapunov exponents in- 
crease as Eq decreases. The increased part of the Lya- 
punov exponents from those at i?o — > oo limit character- 
izes the amplification of the chaotic instability due to the 
interaction between the kicked pendulum and the bulk. 
With this in mind, we define the residual Lyapunov spec- 
tra by 

A\{Eq) = X,{Eq) - lim X,{Eq). (5) 

^From detailed numerical experiments, we find that the 
two scaling relationships; 

Xn-2{Eq) = AXn-2{Eq) cx exp(-/3i/^), (6) 
Ah{Eo) EE ^ AA, cx exp{-f32^/E^) (7) 

i=l 

hold for sufficiently large values of Eq . Here f3i is smaller 
than P2, as shown in Fig.||(a). 

Since A7V-2 = in the limit Eq 00, the scaling of 
Xn-2 represents the decrease in the degree of chaos of 
the kicked pendulum. The sum of the remaining positive 
Lyapunov exponents, then, is expected to correspond to 
the Kolmogorov- Sinai (KS) entropy of the bulk, and the 
increase in the degree of chaos of the bulk due to the 
interaction with the kicked pendulum is given by Ah. 
(Note that the increase of Lyapunov exponents is not due 
to the increase of the average velocity in the bulk. Indeed, 
we have numerically confirmed that the temperature of 
the bulk for finite Eq is not shifted from that in the limit 

Eq 00.) 

Since Ah is associated with the dynamics of the bulk 
and how it is affected by the energy absorbed from the 
kicked pendulum, we expect that it is related to the relax- 
ation time (ifl) Indeed a in Eq.d) and/32 in Eq.(0) 
appears to be equal from our numerical results. The re- 
lationship between {tj^} and Ah is plotted in Fig.^(b), 
which supports the relationship 

Ah cx {tR)-\ (8) 

Thus, we find that the degree of chaos is closely related 
to the relaxation time scale. 

The relationship between the relaxation and the resid- 
ual KS entropy Ah can be roughly explained as follows: 
In a Hamiltonian system, for every orbit, there exists a 
time-reversed orbit, obtained by changing the sign of all 
Pi. However, because of the interaction with the kicked 
element, to obtain the time-reversed orbit of the bulk, 
a slight modification of each pendulum, in addition to 
the simple reversed of each of their momenta, pi — > —pi, 
is needed. Since the KS entropy h is a measure of the 
time rate of loss of the information concerning the initial 
conditions as an orbit evolves, the measure of orbits that 
can be considered the reversed orbit of some given orbit 
over a time t decreases with exp{—ht). Therefore, the 
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time scale for the irreversible relaxation is proportional 
to the inverse of the KS entropy difference between the 
original and the reversed orbits. Hence, the relationship 
Eq.(||) between the time scale for the relaxation and the 
residual KS entropy is reasonable. 

Next, we study dynamics on the long time scale of 
the order of {tn). To compare the relaxation from vari- 
ous kicked energies Eq, it is convenient to use the scaled 
time T = t/{tR), where (tji) is a function of Eq. We note 
that the time scale of relaxation decreases as the energy 
of the kicked pendulum, Kr, decreases. Then the above 
definition of r implies r = t/tii{Kr) along the relaxation, 
where Kr is the coarse-grained energy of the kicked pen- 
dulum over a long time scale and Ir oc exp(a\/ K^)- 

Now, we derive an equation describing the relaxation 
of Kr on the time scale of r. We have computed the dif- 
ference SKr over the time steps St for many samples of 
orbits starting from a given Kr- As shown in Fig.^, the 
data are satisfactory fit by ^^^f^ — —C, with a constant 
C (~ 0.035) for various values of K^ and E. Then the 
coarse-grained equation for the relaxation is obtained as 
dKr/dr = — C. Noting that the interaction energy be- 
tween the bulk and the kicked pendulum is tiny and re- 
calling the conservation of energy, we get dE/dr — C. 
Thus, with regard to the relaxation process, the differ- 
ence in the details of the dynamical properties of the 
system under different conditions are eliminated by con- 
sidering the scaled time r, while the dynamics and prop- 
erties of the bulk seem to be strongly dependent on the 
value of E. ^From the above equation, for the original 
time variable t, Kr decays in the logarithmic time scale. 

Now it is possible to consider the entropy increase of 
the bulk occurring during the relaxation. Since we have 
assumed that the slow change of the bulk part can be de- 
scribed by the slow change of thermodynamic quantities, 
the increase of the entropy of the bulk is estimated as 

dS _dS dE _C 

'd^ ~ dSl^ ~ T' ^ ^ 

where T is the temperature of the bulk. This rela- 
tion is obtained by noting that the process here is suffi- 
ciently slow to define these thermodynamic quantities on 
a coarse-grained time scale. 

Since r = t/tn (x Ah i, by Eq.(||), the relation 

^(xAh (10) 
dt 

for the relaxation process is obtained from Eq.(|9|). The 
proportionality in Eq. (|l^) relates two kinds of directed 
motion towards equilibrium: Ah corresponds to the am- 
plification of chaos in the equi-energy surface of the bulk, 
whereas dS/dt represents the growth of this equi-energy 
surface towards equilibrium. It should also be noted that 
this relationship holds not near the equilibrium state of 
the total system but for the nonlinear relaxation process 
from a highly excited state. 



The relationship between chaotic dynamics and irre- 
versibility has been extensively studied [p"3|-[T5|. A re- 
lationship between the thermodynamic entropy and the 
irreversible information loss was proposed by Sasa and 
Komatsu in the case that an external operation causes 
a transition from one equilibrium state to another. We 
believe that our relation between the residual KS entropy 
and the thermodynamic entropy is related to it, although 
ours applies to the spontaneous relaxation. 

To sum up, we have obtained the following relationship 
between the relaxation time, tn, and the energy supplied 
by an external kick, Eq: t/j cx exp{a\/EQ). Then, this 
relaxation time tn is found to be proportional to the in- 
verse of the residual KS entropy, that is, the difference 
between the KS entropy of the bulk interacting with the 
kicked element and that of the isolated bulk. By rescal- 
ing time by tn, the excited energy is found to relax at 
a rate that is independent of the energy of the kicked 
element and the bulk. Finally, we find evidence that the 
rate of entropy generation is proportional to the residual 
KS entropy. 

Although these three relationships were found through 
numerical study of the Hamiltonian (^, we expect that 
they hold generally for relaxation processes in Hamilto- 
nian systems under the condition that in the limit of high 
excitation, the interaction between the excited elements 
and the unexcited bulk vanishes. Indeed, we have nu- 
merically found that Eqs.(^ and (||) are valid also for a 
coupled pendulum model with the nearest-neighbor cou- 
pling on a square lattice. It is important to examine the 
universality of the three relationships we have found, and 
also to study their relevance in regard to intra-molecular 
relaxation processes, including application to the energy 
transduction carried out by biological molecular motors. 
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FIG. 1. Time series of the momenta Pi {1 < i < N) in 
the process of energy relaxation after an instantaneous kick 
at t = 0. Here, E/N = 0.0205, N = 10 and Eq = 0.35. A is 
the time series of the kicked pendulum. 



FIG. 2. Ensemble average of the relaxation time 
(tR) as a function of the square root of the injected en- 
ergy Eq. Here, (A) E/N = 0.0205 {T = 0.019), (•) 
E/N = 0.0118 (T = 0.008). The dotted line represents 
g{Eo) = exp(15\/2-Eo). In the inset, in order to see the close 
up temperature dependence of (iij), (tR) / g{Eo) is plotted as 
a function of T for four values of Eq: Eq — 0.35 (-I-), 0.4 
(o), 0.45 (x) and 0.5 (□). The data are fit by the function 
13exp(-135T) for smaU T. 



FIG. 3. (a) Dependence of Ajv-2 (o) and Ah (□) on Eq 
for E/N = 0.0205 (T = 0.019). The data are fit by the 
lines AAjv-2 oc exp(4.1\/2£'o) and Aft oc exp(15\/2-£'o). (b) 
{tR) as a function of Aft for E/N = 0.0205 (T = 0.019). 
To calculate the Lyapunov Spectrum Ai(_Eo), we adopt the 
criterion \Eo - Kr\ < S, with S = 0.00375. The calculated 
Ai(_Eo) is found to converge for small S. 

FIG. 4. {5Kr)/5r for 5t = 0.01 (o) and 0.001 (□). upper: 
dependence on the temperature T. lower: dependence on the 
energy Kr of the kicked pendulum. The data suggest that 
{5Kr)/ST is independent of St, Kr and T. 
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